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ABSTRACT 

Context. Solar prominences are partially ionised plasmas displaying flows and oscilla- 
tions. These oscillations show time and spatial damping and, commonly, have been 
explained in terms of magnetohydrodynamic (MHD) waves. 

Aims. We study the spatial damping of linear non-adiabatic MHD waves in a flowing 
partially ionised plasma, having prominence-like physical properties. 
Methods. We consider single fluid equations for a partially ionised hydrogen plasma 
including in the energy equation optically thin radiation, thermal conduction by elec- 
trons and neutrals, and heating. Keeping uj real and fixed, we have solved the obtained 
dispersion relations for the complex wavenumber, k, and have analysed the behaviour 
of the damping length, wavelength and the ratio of the damping length to the wave- 
length, versus period, for Alfven, fast, slow and thermal waves. 

Results. The results indicate that in presence of a background flow new strongly damped 
fast and Alfven waves appear which depend on the joint action of flow and resistivity. 
The damping lengths of adiabatic fast and slow waves are strongly affected by partial 
ionisation, which also modifies the ratio between damping lengths and wavelengths. 
Moreover, the behaviour of adiabatic fast waves looks very similar to that of Alfven 
waves. For non-adiabatic slow waves, the unfolding in wavelength and damping length 
induced by the flow, allows that an efficient damping can be found for periods compat- 
ible with those observed in prominence oscillations. This effect is enhanced when low 
ionised plasmas are considered. 

Conclusions. Since flows are ubiquitous in prominences, in the case of non-adiabatic 
slow waves and within the range of periods of interest for prominence oscillations, the 
joint effect of both flow and partial ionisation leads to a ratio between damping length 
and wavelength denoting a very efficient spatial damping. For fast and Alfven waves, 
the most efficient damping occurs at very short periods not compatible with those 
observed in prominence oscillations. 
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1. Introduction 

A solar prominence is a cool (T ~ 10* K) and dense (p ~ 5 • 10^^ kg/m"^) mass of gas 
located in the much denser and hotter solar corona. Although it is not well understood 
how this structure can live for a long time within the solar corona, it is thought that its 
support and thermal isolation are of magnetic origin. Thanks to ground and space-based 
observations, evidence about the presence of small-amplitude oscillations in prominences 
and filaments is now widely available and detailed information can be found in Oliver and 
Ballester (2002), Banerjee et al. (2007), Oliver, (2009), and Mackay et al. (2009). The os- 
cillations are mainly detected through the periodic Doppler shifts of spectral lines, with 
periods from a few minutes to hours, and the observations have shown that they are of 
local nature, that they can be coherent over large regions of prominences/filaments, that 
simultaneously flowing and oscillating features are present, and that oscillations display 
time and spatial damping (Terradas et al. 2002). These small-amplitude oscillations have 
been commonly interpreted in terms of linear magnetohydrodynamic (MHD) waves and 
the damping of oscillations has been studied by considering different dissipative mecha- 
nisms. For instance, in the case of fully ionised plasmas the time damping of prominence 
oscillations has been studied by considering non-adiabatic MHD waves, including the effect 
of optically thin radiation, thermal conduction and heating, in unbounded and bounded 
prominence plasmas (Carbonell et al. 2004; Terradas et al. 2005); furthermore, assuming 
that the threads that seem to compose solar filaments can be modeled as cylindrical flux 
tubes, the time damping of homogeneous and inhomogeneous cylindrical flux tubes, with 
prominence physical conditions, has been modeled using thermal mechanisms and resonant 
absorption, respectively (Soler et al. 2008; Arregui et al. 2008; Soler et al. 2009b). 

A typical feature in prominence oscillations is the presence of flows which are observed 
in Ha, UV and EUV lines (Labrosse et al. 2009). In Ha quiescent filaments, the observed 
velocities range from 5 to 20 km/s (Zirker et al. 1998; Lin et al. 2003, 2007) and, because 
of physical conditions in filament plasma, they seem to be field-aligned. In the case of 
active region prominences, fiow speeds can be higher. Recently, observations made with 
Hinode/SOT by Okamoto et al. (2007) reported the presence of synchronous vertical oscil- 
latory motions in the threads of an active region prominence, together with the presence 
of flows along the same threads. However, in limb prominences diff'erent kinds of flows are 
observed and, for instance, observations made by Berger et al. (2008) with Hinode/SOT 
have revealed a complex dynamics with vertical downflows and upflows. Taking all this into 
account, Carbonell et al (2009) explored the time damping of non-adiabatic slow and ther- 
mal waves in an unbounded prominence medium with a background flow, while Soler et al. 
(2008, 2009a) investigated the time damping of the oscillations of an individual prominence 
thread and of a threaded prominence when both mass flows and non-adiabatic processes 
are considered. 
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Although there is a large amount of observational evidence about the time damping of 
MHD waves propagating in coronal structures such as coronal loops or prominences, current 
observational information about the spatial damping of MHD waves in coronal structures 
is still scarce. Time damping is usually produced when an impulsive perturbation excites a 

medium and waves start to propagate, then, the wave amplitude decreases as the time goes 
by. Alternatively, if the medium is excited by a continuous driver with a fixed frequency, the 
spatial damping is detected by the decrease of the amplitude when the wave propagates. 
Since, probably, most of the phenomena exciting waves in solar structures are of impulsive 
origin rather than produced by continuous drivers, time damping of waves is more often 
observed. 

Concerning spatial damping, Dc Moortcl ct al. (2002) and Dc Moortcl & Hood (2003, 
2004) , using thermal conduction, compressive viscosity, gravitational stratification and field 
line divergence, studied the spatial damping of driven and non-driven slow MHD waves in 
coronal conditions, applying the obtained results to the case of standing and propagating 
slow waves in coronal loops observed with SOHO and TRACE. In the case of prominences, 
Terradas et al. (2002) analyzed small amplitude oscillations in a polar crown prominence 
reporting the presence of a plane propagating wave as well as an standing wave. The plane 
wave propagates in opposite directions with wavelengths of 67, 500 and 50, 000 km and 
phase speeds of 15 km s~^ and 12 km s~^, while in the case of the standing wave, the 
estimated wavelength is of 44,000 km and the phase speed of 12 km s~^. These authors 
also reported that in the case of the propagating wave, which was interpreted as a slow 
MHD wave, the amplitude of the oscillations spatially decreases in a substantial way in a 
distance of 2 — 5 x 10"* km from the location where wave motion was being generated. This 
distance can be considered as a typical spatial damping length, L^, of the oscillations. 

Prom the theoretical point of view, and assuming a fully ionised plasma, Ballai (2003) 
performed a qualitative study of the spatial damping of linear compressional waves in 
solar prominences by considering different dissipative mechanisms such as isotropic and 
anisotropic viscosity, isotropic magnetic diffusivity, isotropic radiative damping (Newton's 
law) and anisotropic thermal conduction. The conclusion was that thermal radiation can 
damp compressional waves as well as that waves can be also damped by anisotropic thermal 
conduction provided that they have enough short wavelengths. Using linear non-adiabatic 
MHD waves, which included optically thin radiation, thermal conduction and heating, 
Carbonell et al. (2006) made a quantitative study of the spatial damping of slow, fast 
and thermal waves in a fully ionised unbounded prominence medium. In the frequency 
space, they determined the regions where radiation or thermal conduction are the dominant 
damping mechanisms, the critical frequencies at which the dominance changes from one 
mechanism to another, and the fact that different heating mechanisms do not strongly affect 
the damping. The most important conclusions were that the thermal wave is propagating 
but strongly damped, and that slow waves are efficiently damped by thermal effects, which 
docs not happen for fast waves. Singh et al. (2006) did a similar study but considering only 
Newtonian radiation and, later, Singh et al. (2007) considered again the same problem 
including Newtonian radiation and turbulent viscosity. 
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The typical composition of solar prominences is 90% hydrogen and 10% helium and they 
constitute partially ionised plasmas since hydrogen lines are observed (Labrosse et al. 2009). 
The exact ionisation degree of prominences is unknown and the ratio of electron density to 
neutral hydrogen density seems to be in the interval 0.1 — 10 (Patsourakos and Vial, 2002). 
From a theoretical point of view, partial ionisation was already considered by Mcrcicr and 
Heyvaerts (1977) when they studied the difussion of neutral atoms due to gravity; also, 
Gilbert ct al. (2002) studied the diffusion of neutral atoms in a partially ionized prominence 
plasma, concluding that the loss time scale is much longer for hydrogen than for helium. 
Recently, Gilbert et al. (2007) have investigated the temporal and spatial variations of the 
relative abundance of helium with respect to hydrogen in a sample of filaments. They have 
found that a majority of filaments show a deficit of helium in the top part while in the 
bottom part there is an excess. This seems to be due to the shorter loss time scale of neutral 
helium compared to that of neutral hydrogen. The consideration of prominences as made 
of partially ionised plasma is extremely important for the physics of prominences, and the 
effects on MHD waves in prominences need to be taken into account. From the theoretical 
point of view, and in the framework of laboratory plasma physics, a extense literature 
about wave propagation in a partially ionized multifluid plasma can be found (Watanabe, 
1961a, b; Tanenbaum, 1961; Tancnbaum and Mintzcr, 1962; Woods, 1962; Kulsrud and 
Pcarcc. 1969; Watts and Hanna, 2004). In astrophysical plasmas, the typical frequency of 
MHD waves is much smaller than the coUisional frequencies between species. In such a 
case the single fluid approach is usually adopted and it has been applied to wave damping 
in the solar atmosphere (De Pontieu, et al. 2001; Khodachenko et al. 2004; Leake et al. 
2005). In the case of solar prominences, Forteza et al. (2007) derived the full set of MHD 
equations for a partially ionized, single fluid plasma and applied them to study the time 
damping of linear, adiabatic waves in an unbounded prominence medium. Later, this study 
was extended to the non-adiabatic case by including thermal conduction by neutrals and 
electrons, radiative losses and heating (Forteza et al. 2008). Because of the effect of neutrals, 
in particular that of ion-neutral collisions, a generalized Ohm's law has to be considered, 
which causes some additional terms to appear in the resistive magnetic induction equation 
in comparison to the fully ionized case. Among these additional terms, the dominant one in 
the linear regime is the so-called ambipolar magnetic diffusion, which enhances magnetic 
diffusion across magnetic field lines. Furthermore, one of the interesting effects produced 
by the consideration of partial ionisation and ion-neutral collisions is that Alfven waves 
can be damped, which cannot be obtained by means of thermal mechanisms, and Singh & 
Krishnan (2009) have studied the time damping of Alfven-like waves in a partially ionised 
plasma representing a particular model of the solar atmosphere. For bounded media, Soler 
et al. (2009c) have applied this formalism to the study of the time damping of fast, Alfven 
and slow waves in a partially ionised filament thread modeled as a cylinder, and Soler et 
al. (2009d) have used a cylindrical filament thread, having an inhomogeneous transition 
layer between prominence and coronal material, to study the influence of partial ionisation 
on the time damping of fast kink waves which is caused by resonant absorption in the 
inhomogeneous layer. 
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Up to now, and to our knowledge, no study of the spatial damping of MHD waves 
in a partially ionised plasma has been performed. Therefore, since solar prominences are 
partially ionised plasmas in which material flows and oscillations are present and since 
commonly these oscillations are interpreted in terms of MHD waves, our main aim here 
is to explore the theoretical and observational effects associated to the spatial damping of 
non-adiabatic MHD waves in an unbounded and partially ionised plasma, with prominence- 
like physical conditions, when a background flow is present. The layout of the paper is as 
follows: In Sect. 2, the equilibrium model and some theoretical considerations are presented; 
in Sect. 3, the spatial damping of Alfven waves is studied; in Sect. 4, the spatial damping 
of magnetoacoustic waves is considered; finally, in Sect. 5, conclusions are drawn. 

2. Model and Methods 

2.1. Equilibrium Properties 

As a background model, we use a homogeneous unbounded medium threaded by a uni- 
form magnetic field along the a;-direction, and with a field-aligned background flow. The 
equilibrium magnitudes of the medium are given by 



with Bq and vq constants, and the effect of gravity has been ignored. Since we consider 
a medium with physical properties akin to those of a solar prominence, the density is 
po = 5 X 10-11 kg/m^, the temperature Tq = 8000 K, the magnetic fleld |Bo|= 10 G, 
and, in general, a fleld-aligned flow with vq = 10 km/s, simulating the typical flows found 
in quiescent filaments, has been considered. 

2.2. Basic and Linearised Equations 

The derivation of the single-fluid MHD equations for a partially ionised hydrogen plasma 
can be found in Forteza et al. (2007). Later, these equations were modified in order to in- 
clude non-adiabatic processes in the energy equation (Forteza et al. 2008), and the physical 
meaning of all the terms and quantities \ised in the following can be found in those papers. 
When a background flow is considered, the single fluid basic equations for the study of 
non-adiabatic MHD waves in a partially ionised plasma can be derived from Forteza et al. 
(2007, 2008) and Carbonell et al. (2009). The most important difference with respect to 
the non-adiabatic case without a background flow is that instead of the operator ^ we will 
have ^ -|- vo • V in the MHD equations (Goedbloed and Poedts, 2004). Before proceeding, 
some remarks are in order: We consider a partially ionised hydrogen plasma characterized 
by a plasma density, po, temperature, Tq, and number densities of neutrals, n„, ions, rij, 
and electrons Ug, with rig = n,. Thus, the gas pressure is po = (2n, -|- n„)fcBTo, where ks 
is Boltzmann's constant. The relative densities of neutrals, and ions, ^j, are given by 



Po — const., Po = const., Tq = const.. 



Bo = Box, vo = Vox, 



'71 



6 = 



(1) 



rii -\-n- 



"n 



"n 
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where we have neglected the contribution of electrons. We can now define a ionization 
fraction, /t, which gives us information about the plasma degree of ionization, 

* = T^- 

and the quantity 

where a„ is a friction coefficient (Braginskii, 1965; Khodachenko et al. 2004; Leake et al. 
2005). For a fully ionised plasma jl = 0.5 while for a neutral plasma /t = 1. Depending 
on the value given to /i and to Spitzcr's (rj) and Cowling's (rjc) resistivities, we may have 
different types of plasmas: A fully ionised ideal plasma (FIIP) is characterised by /2 = 0.5, 
Tj = rjc = ^ — 0; a fully ionised resistive plasma (FIRP) is characterised by /x = 0.5, t] = r]c 
and S = 0; and a partially ionised plasma (PIP) is characterised by 0.5 < /i < 1, r; ^ r]c 
and S ^ 0. Although our main aim is the study of spatial damping of MHD waves in PIP, 
the case of FIRP will be also considered since it has received almost no attention in the 
literature. Another important issue for our study is the numerical value and behaviour of 
the sound (cg) and Alfven speeds {vg). For fixed the density and magnetic field, the Alfven 
speed has a constant numerical value of 126.15 km/s, however, since the sound speed 
depends on gas pressure, which is a function of the number densities of ions and neutrals, 
its numerical value is not constant but depends on the ionisation fraction considered. For 
a fully ionised plasma, the sound speed is 14.84 km/s while for a partially ionised plasma 
with fl = 0.95 its value decreases to 10.76 km/s. This variation in the sound speed can be 
important since depending on the flow speed and ionisation fraction chosen, the flow speed 
could be greater than, smaller than or equal to the sound speed which affects the direction 
of propagation of slow and thermal waves (Carbonell et al. 2009). In our case, the flow 
speed is subsonic and subalfvenic which seems to be a typical feature of flows observed in 
quiescent fllaments. 

Next, to obtain the dispersion relation for linear MHD waves, we consider small per- 
turbations from the equilibrium of the form 

B{t,r) = Bo + Bi{t,r), p{t,r) = po + Pi(t,r), 
p{t, r) = po + pi{t, r), T{t, r) = To + T^it, r), 
v{t,r) = vo +vi{t,r). 

and we linearise the single fluid basic equations. Since the medium is unbounded, we 
perform a Fourier analysis in terms of plane waves and assume that perturbations behave 
as 

/i(r,t) = /e*(-*-'=-^). (4) 

and with no loss of generality, we choose the z-axis so that the wavevector k lies in the 
xz-planc (k = k^x + k^z), and introducing the propagation angle, 9, between k and Bq, 
the wavcnumbcT components can be expressed as k^ — kcos6 and k^ — ksinO. With this 
approach, the operator ^ + vq • V becomes i{u — k^vo), which points out that in the 
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presence of a background flow the frequency has a Doppler shift given by kxVo and that 
the wave frequency, w, for the non-adiabatic case with flow can be obtained from 

w = fi + kxVQ, (5) 

Q being the wave frequency for the non-adiabatic case without flow. Also, these frequencies 
can be described in a different manner, fl corresponds to the frequency measured by an 
observer hnked to the flow inertial rest frame, while uj would correspond to the frequency 
measured by an observer linked to an external inertial rest frame. 

After the Fourier analysis, the following linearised scalar equations are obtained. 

^Pi - PoikxVix + kzViz) = 0, (6) 

^Povix - kxPi = 0, (7) 

^Poviy + ^^kxBiy = 0, (8) 

ftpovi^ - k^pi -] -{kxBi^ - k^Bix) ^0, (9) 

M 

iO. {pi - CgPi) + (7 - 1) (Ko||fc^ + Knfc^ + PoLt) Ti + 

+(7-l)(L + poip)pi =0, (10) 

Bix {iil + kit] + kf rjc) + (r) ~ Vc)kxkzBiz - 

-Boxkz{iviz~ k^'E.pi) ^Q, (11) 

Biy {in + kli]c + klrj) + iBoxkxViy = 0, (12) 

Biz {i^i + k^rjc + kit]) + [rj - r]c)kxkzBix + 

+Boxkx{iviz-kz'E.pi)^Q, (13) 

kxBix + kzBiz = 0, (14) 

El-Pl-Il^O. (15) 
Po Po To 

2.3. Dispersion relation for Alfven waves 

Equations ([5]) and (|12p are decoupled from the rest and from them we can obtain a disper- 
sion relation for Alfven waves in a partially ionised plasma with a background flow, which 
is, 

n'^ - iQk'^irjc cos^ e + rism^ 9) - vlk^ cos^ 9 = (16) 

where 9 is the angle between the wavenumber vector and the magnetic field. In absence of 
flow, Eq. (fT6)) reduces to, 

r(6l)2 C0S2 

where T{9) is a modified and complex Alfven speed {T{9) = Fj- + iFi) defined as 

T{9f = vl+ iuj (77c + V tan^ &) ■ (18) 

which was introduced by Forteza et al. (2008). In the case of a fully ionised ideal plasma 
without a background flow Eq. (jl6p reduces to the well known expression for Alfven waves. 
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2.4. Dispersion relation for magnetoacoustic waves 

From the remaining linearised equations, imposing that the determinant of the algebraic 
system be zero, we obtain our general dispersion relation for thermal and magnetoacoustic 
waves in presence of a background flow, which is given by, 

{n^ ~ k^A^){ierjcn - n^) + k^Hn" - klA^) + 

+ik^klvlA^Epon = (19) 

where is the non-adiabatic sound speed squared (Forteza et al. 2008, Soler et al. 2008) 
and is defined as 
^ ^A-H + icln 

= ^ (20) 

pa 

where A and H, including optically thin radiative losses, thermal conduction by electrons 
and neutrals, and a constant heating per unit volume, are defined in Forteza et al. (2008). 
When Eq. (|19[) is expanded, it becomes a seventh degree polynomial in the wavenumber k. 
By setting the appropiate vakies to the corresponding quantities, the dispersion relation 
l|19p is consistent with the case of an adiabatic and partially ionised plasma without flow 
(Forteza et al. 2007), the case of a non-adiabatic and partially ionised plasma without 
flow (Forteza et al. 2008) and the case of a non-adiabatic, fully ionised plasma with a 
background flow (Carbonell et al. 2009). 

3. Spatial damping of Alfven waves in a partially ionised plasma 

Since we are interested in the spatial damping of MHD waves we consider the frequency, 
uj, to be real and seek for complex solutions of the wavenumber k expressed as fc = fcr + iki. 
The wavelength of the waves is given by A = j^, the damping length by Ld = p and 



the damping length per wavelength is and, in general, we consider a propagation angle 



3.1. Spatial damping of Alfven waves without background flow 

In this case, our governing dispersion relation is given by Eq. (jl7p . and the wavenumbers 
are 



k=±,l , , (21) 

y w2 -t-ia;(?7c + r?tan2 6') 

representing two Alfven waves propagating in opposite directions. The real part (fc^) of 
these wavenumbers is, 

kr^oj sec 9 ^ (22) 

while the imaginary part (fc;) is, 

— sec 6 
^2 [w4 4-a;2(,yc +7?tan2 9)] 
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Fig. 1. Left panels: Damping length, wavelength, and ratio of the damping length to the 
wavelength versus the period for Alfven waves in a FIRP (solid) and in PIP with jl = 0.8 
(dashed), fi = 0.95 (dotted), and p, = 0.99 (dash-dotted). Right panels: Damping length, 
wavelength, and ratio of the damping length to the wavelength versus period for the three 
(solid, dashed, dotted) Alfven waves in a PIP with jj, = 0.8 and with a background flow 
of 10 km/s. In all the Figures, the shaded region corresponds to the interval of observed 
periods in prominence oscillations. 



(23) 



Figure [T] (Left panels) shows a plot of the damping length, wavelength, and the ratio of 
the damping length to wavelength, versus the period, for Alfven waves. The plots have 
been made for four different ionisation fractions and the shaded region corresponds to the 
interval of observed periods in prominence oscillations. When a period within the shaded 
region is considered, we observe that the damping length decreases in a substantial way 
when the amount of neutrals in the plasma increases. Then, when ion-neutral collisions 
are present the spatial damping of Alfven waves is enhanced for periods greater than 1 s. 
Also, we can observe that for a FIRP the behaviour of the damping length, wavelength 
and their ratio versus period is linear, except for periods below 10^^ s. However, when a 
PIP is considered a deviation from the linear behaviour appears for periods below 1 s. This 
is due to the joint effect of the terms including frequency and resistivities in the real (kr) 
and imaginary (fc^) parts of the wavenumber k since going to shorter periods frequency 
increases and its role becomes more important. The spatial damping of Alfven waves in 
PIP with physical properties like those of quiescent prominences is very efficient for periods 
below 1 s, however, within the interval of periods of interest for prominence oscillations, it 
is only efficient, (-^ ~ 1 — 10), when almost neutral plasmas are considered. 
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3.2. Spatial damping of Alfven waves with background flow 

Now, our dispersion relation is given by Eq. (jl6p which once expanded becomes a cubic 
polynomial in the wavenumber fc, such as, 

Wo ('7c cos^ 6 + ?7sin^ 6)k^ + 
~^ [(^0 ~ "^a ~ ivc^) COS 9 — irjuj sin 9 tan 0] fc^ — 2ujvok + 

sec 61 = 0, (24) 

The increase in the degree of the dispersion relation, with respect to the case without flow, 
is produced by the joint presence of flow and resistivities. In this case, we should obtain 
three propagating Alfven waves, therefore. Figure [T] (Right panels) shows the numerical 
solution of dispersion relation ([M)) for the three Alfven waves in a partially ionised plasma 
with a background flow. For all the interval of periods considered, a strongly damped third 
Alfven wave appears, while on the contrary, and like in 13. 1[ the other two Alfven waves 
are very efficiently damped for periods below 1 s. However, within the interval of periods 
typically observed in prominence oscillations these waves are only efficiently attenuated 
when almost neutral plasmas are considered. Furthermore, the following approximations 
for the different wavenumbers corresponding to two of the expected Alfven waves can be 
calculated. 



Y (wo ± Wa)2 + za;(77c + ?7tan2 6*) 
whose real part (fcr) is given by. 



J {vq ± v^)^ + y^vo ± Wa)^ + uj^ {vc + tau^ 9) 

uj sec 9^ ^ (26) 

^2 [{vo ± Va)4 + ^Hvc + '/tan2 9)] 



while the imaginary part (fc;) is, 

— in'^ spp ft 



2[(z;o±i;a)4+w2(r/c+r7tan2 9)] 

ivc + 1] tan^ 9) 
(wo ± Wa)2 + v/(^o ± ^^a)^ + ^"^{ric + ?? tan^ 



(27) 



From the above expressions, if we consider a FIIP we recover the dispersion relation for 
Alfven waves with a background flow (Carbonell et al. 2009), and if we remove the flow, 
the well-known dispersion relation for Alfven waves is recovered. Since the flow speed is 
much smaller than the Alfven speed, the effect of the flow on the real and imaginary parts 
of the wavenumber is very small, then, the wavelengths and damping lengths are similar 
to those in 13.11 The third remaining wavenumber of Eq. (j24p can be approximated by, 

fc = , . (28 

Wo cos ft va {f]c cos^ ft + i] sin^ ft) 

corresponding to the third Alfven wave. All the above analytical approximations display 
an excellent agreement with the numerical results, and the presence of the third Alfven 
wave, given by Eq. (j28p . fully depends on the join presence of flow and resistivities since. 
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otherwise, the dispersion relation ([M)) would be quadratic. For an external observer to 
the flowing plasma, this additional wave could be detectable, although its strong spatial 
damping would make its detection very difficult. For an observer linked to the flow inertial 
rest frame, only the two usual Alfven waves, modified by resistivities, would be detected. 

4. Spatial damping of Magnetoacoustic waves in a partially ionised plasma 

Our general dispersion relation for non-adiabatic magnetoacoustic waves in presence of a 
background flow is given by Eq. (jl9p . Because of the complexity of this dispersion relation, 
and in order to be able to understand the results, we have split our study in a sequence 
of four different cases having dispersion relations of increasing complexity. In the flrst two 
cases, the behaviour of adiabatic magnetoacoustic waves in a non-flowing (|4.ip and flowing 
plasma (|4.2p is considered; while in the last two cases, the behaviour of non-adiabatic waves 
in a non-flowing (14. 3p and flowing plasma (j4.4p is studied. 

4.1. Adiabatic Magnetoacoustic waves without background flow 

Setting A = H — in Eq. (l20l) . wo = in expression ([5]) and substituting in Eq. (|T9l) . 
we obtain the dispersion relation for adiabatic magnetoacoustic waves in a PIP without a 
background flow, which is 

(w^ - k'^cl){ik^r]cLo - u?) + k^vl{uj'^ - klcl) + ik^ klvlc^Epouj = (29) 
already found by Forteza et al. (2007) 

4.1.1. Fully ionised resistive plasma 

Now, imposing FIRP conditions from 12. 2[ the dispersion relation simplifles to, 

(w^ - ecDiik^riu - ^2) + k^vliu'^ - klcl) = (30) 

which once expanded gives a fourth degree polynomial in the wavenumber k, 

fc*c^(cos^ 0vl + iujri) ~ k^uj^icl + vl + iujr]) -|- w"* = 0. (31) 

Furthermore, when only longitudinal propagation is allowed {6 = 0), dispersion rela- 
tion can be factorized as 

(u;2 _ k^d) [u^ - k^ivl + z^,/)] = (32) 
giving place to two undamped slow waves with a dispersion relation given by, 

k^=4 

c 

and, since for longitudinal propagation fast waves become Alfven waves, we obtain two 
Alfven waves, damped by resistivity, whose dispersion relation is. 



For these Alfven waves, the real (fc^) and imaginary (ki) parts of the wavenumber k can be 
obtained from expressions and setting = 0. Then, in the case of a FIRP, when 
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oblique propagation is not allowed, resistivity does not affect slow vifaves, which propagate 
undamped, and only affects fast waves. Next, if we allow oblique propagation {0 < 9 < tt/2) 
in dispersion relation (j3ip , the solutions for the wavenumbers are given by, 

.2 ^^ivl+(i+iV^) , 



2iCs{r]uj — iv^ cos^ 



■^V(^a + + iy^y - 4icg(?7cj - iv l cos2 6) ^^^^ 



and 



2iCg(r]ijj — iv"^ cos^ 9) 
,2 + 



2iCg{r]uj — iv^ cos^ 9) 



^V('"a + eg + »??^)^ - 4icg(?7cj - ivl cos^ 9) 
2iCg{riu! — ivi cos^ 9) 

where fc^ and fc^ are the wavenumbers corresponding to two coupled and damped fast and 
slow waves propagating in opposite directions. Figure [H displays the behaviour of damping 
length, wavelength and the ratio of damping length to wavelength versus period for fast 
and slow waves, respectively. In these figures it is clearly observed that the behaviour of 
the wavelength is the same for a FIRP and a FIIP, and that the only difference occurs at 
periods below 10~^ s due to resistivity. However, looking to the ratio of the damping length 
to the wavelength, we can conclude that in a FIRP, and within the interval of periods of 
interest, the spatial damping of both waves is negligible. 



4.1.2. Partially ionised plasma 

The dispersion relation is given by Eq. (I29p , and considering only longitudinal propagation 
we obtain, 

(w^ - k^cl) [a;2 _ fc2(^2 ^ ^ Q (35) 

This could suggest that the consideration of longitudinal propagation in a FIRP or in a 
PIP leads to the same dispersion relation, however, there is an important difference, while 
for FIRP both resistivities have the same numerical value, for PIP the numerical value 
of Cowling's resistivity is much greater than that of Spitzer's resistivity. Furthermore, for 
longitudinal propagation the slow waves are not influenced by Cowling's resistivity. Once 
expanded, Eq. (|35|) gives the following fourth degree polynomial in the wavenumber fc. 



fc^'iCg — ivl^ cos^ 9 + i;^po'^'^(cos^ ^ — 1)] — 

-k'^uj'^{vl + c^^^'iricuj)+uj'^ = Q. (36) 

After solving this biquadratic dispersion relation, we obtain, 

- ^ (37) 

*2 - ^ (38) 

with 

^ = (^'a + + im^) 

B ^ ivl + cl + wycw)" - 
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Fig. 2. Damping length, wavelength, and ratio of the damping length to the wavelength 
versus period for the adiabatic fast (left panels) and slow (right panels) waves in a FIIP 
(solid), a FIRP (dashed) and a PIP with fl ^ 0.8 (dotted) and fl = 0.95 (dash-dotted). 
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Fig. 3. Damping length, wavelength, and ratio of the damping length to the wavelength 
versus period for the adiabatic fast (left panels) and slow (right panels) waves in a FIIP 
(solid) and in a PIP with /i = 0.8 (dashed). The background flow speed is 10 km/s. 

-4iCg [rjcLJ - ivl cos^ 9 - vlpoE{cos^ 9 - 1)] 
C = 2icl [ricijJ - ivl cos^ 9 - w^poS(cos^ 9 - 1)] 



and the wavenumbers are, 



ki = ±Vw2fc2^ 

fc2 =±V^ 



(39) 
(40) 
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versus period for the non-adiabatic fast (left panels), slow (right panels) waves in a FIIP 
(solid), a FIRP (dashed), and a PIP with jl = 0.8 (dotted) and jl = 0.95 (dash-dotted). 



Figure [2] displays the behaviour of damping length, wavelength and the ratio of damping 
length to wavelength versus period for fast and slow waves in a PIP, respectively. The 
damping length of both slow and fast waves is severely influenced by ion-neutral collisions 
showing a strong dependence on the period for periods greater than 1 s while for shorter 
periods the dependence becomes weaker. This figure also shows that with respect to a FIIP 
the wavelength of the fast waves is slightly affected by the partial ionisation, deviating from 
the linear behaviour for periods below 1 s, while the wavelength of slow waves is not affected 
at all. Within the interval of observed periods in prominence oscillations, when ionisation 
decreases the ratio between the damping length and the wavelength also decreases for both 
waves and the spatial damping becomes more efficient. The maximum efficiency of the 
spatial damping for fast waves is attained for periods below 1 s while for slow waves the 
maximum of efficiency is attained at a period which depends on the ionisation fraction. The 
location of this maximum moves towards long periods when the ionisation of the plasma 
decreases, but when almost neutral plasmas are considered it is still located at a period 
slightly greater than 1 s, ouside the region of interest. Finally, comparing Figures [T] and [2] 
it becomes obvious that the behaviour of Alfven and fast waves is quite similar. 

4.2. Adiabatic Magnetoacoustic waves with background flow 

Setting A — H — in Eq. (|20p and substituting in Eq. we obtain the dispersion 

relation for adiabatic magnetoacoustic waves in a PIP with a background flow, which is 



(41) 
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Fig. 5. Damping length, wavelength, and ratio of the damping length to the wavelength 
versus period for the non-adiabatic thermal wave in a FIIP (solid), a FIRP (dashed), and 
a PIP with p. = 0.8 (dotted) and p, = 0.95 (dash-dotted). 

4.2.1. Fully ionised resistive plasma 

Considering FIRP conditions from 1 2. 2 [ the dispersion relation (PT|) becomes, 

{n^ - k^cl){ik^r^VL - k^vKn^ - klcl) = (42) 

The dispersion relation is a fifth degree polynomial in the wavenumber k, and for this reason 
we expect two slow waves and two fast waves which, for the flow speed considered, propagate 
in opposite directions, plus an additional wave. When only longitudinal propagation is 
considered, the above dispersion relation becomes, 

{^f -k^cl)[k'^{i'nVl^vl)-n^]=Q (43) 

and slow waves are decoupled from fast waves propagating undamped, while fast waves are 
damped by resistivity. The wavenumbers corresponding to the undamped slow waves are 
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Fig. 6. Damping length, wavelength, and ratio of the damping length to the wavelength 
versus period for the non-adiabatic fast (left panels), slow (right panels) waves, without 
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Fig. 7. Damping length, wavelength, and ratio of the damping length to the wavelength 
versus period for the non-adiabatic fast (left panels), slow (right panels) waves in a PIP {fl — 
0.8) without neutrals thermal conduction (solid); without electronic thermal conduction 
(dashed) and without thermal conduction (dotted); in a FIIP without thermal conduction 
(dash-dotted), and in a FIRP without thermal conduction (long-dashed). 



given by, 
k " 



, (44) 

Wo ± Cs 

Concerning the fast waves, now Alfven waves because of longitudinal propagation, the 
corresponding dispersion relation, given by the second factor in Eq. (03]), is equivalent 
to Eq. (1^^ when 9 = Q. The solutions to this dispersion relation are given by Eqs. ([^5]) 
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and (j28p with 9 — 0, and we obtain three Alfven waves similar to those studied in 13.21 
Then, the expected additional wave mentioned above is a fast wave which, for longitudinal 
propagation and when a background flow is present, becomes the third Alfven wave already 
found in 13.21 When oblique propagation is allowed, fast and slow waves become coupled 
and the dispersion relation (j42l) is solved numerically. Figure [3] displays the behaviour of 
the damping length, wavelength and the ratio of the damping length to wavelength versus 
period for fast and slow waves in a FIRP. Because of the strong difference between Alfven 
and flow speeds, in the case of fast waves the unfolding in wavelength and damping length 
due to the flow is not evident, while for slow waves it is clearly seen. In both cases the 
behaviour of the wavelength versus period is linear, similar to what happens for a FIIP and 
the only difference is provided by the unfolding produced by the flow. The damping length 
also behaves linearly with period and, for slow waves, the unfolding in wavelength and 
damping length produces two different curves for the ratio L^/X. In both curves, the most 
efficient spatial damping appears for periods far away from those of interest in prominence 
oscillations. For fast waves, the ratio L^/X behaves linearly with period and its value is 
very large within the region of periods of interest. The behaviour of the third fast wave is 
quite different from the other two fast waves and very similar to that of the third Alfven 
wave shown in l3.2[ being strongly damped within the interval of periods considered. 

4.2.2. Partially ionised plasma 

Our dispersion relation is given by Eq. (HT|) which is a fifth degree polynomial in fc, and 
when longitudinal propagation is considered we recover the results of 14.2.11 with slight 
differences due to the different numerical value of Cowling's resistivity. Such as it is shown 
in Figure 131 at periods longer than 0.1 s, the damping lengths of both slow and fast waves 
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increase linearly with the period. However, for periods below 0.1 s the damping length 
slowly decreases in the case of fast waves and becomes constant for slow waves. Furthermore, 
the wavelength of fast waves for periods below 0.1 s increases with respect to the FIRP case 
while the wavelengths corresponding to slow waves are only slightly modified. Concerning 
the ratio of the damping length versus the wavelength, the unfolding in wavelength caused 
by the fiow and the change of the damping length due to the partial ionisation produces 
that fast waves are much more efficiently attenuated, than in a a FIRP, for any period, but 
especially at periods below 0.1 s, while for slow waves the peak of maximum efficiency is 
displaced towards long periods when ionisation decreases, and for almost neutral plasmas 
would approach the region of periods usually observed in prominence oscillations. Again, 
the behaviour of the remaining third fast wave is very similar to that found for the third 
Alfven wave discussed in 13.21 

In absence of flow the dispersion relation (|42|) becomes a fourth degree polynomial in k 
and the third fast wave disappears. Such as it happens for Alfven waves with a background 
flow, this third fast wave appears due to the joint action of flow and resistivity, since when 
a Flip with a background flow is considered, the dispersion relation also becomes a fourth 
degree polynomial in k and the third fast wave is absent. 



4.3. Non-adiabatic Magnetoacoustic waves in a prominence plasma witliout bacl<ground flow 

Setting uq = in expression ([5]) and substituting in Eq. (I19L we obtain the dispersion 
relation for non-adiabatic magnetoacoustic waves in a PIP without a background flow, 
which is 



+ik^klvlA^Epouj = 



(45) 
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4.3.1. Fully ionised resistive plasma 

Once imposed the conditions corresponding to a FIRP, the foUowing dispersion relation, a 
sixth degree polynomial in the wavenumber k, is obtained, 

(uj^ - eA^){ik^7juj - uj^) + k\l{uo^ - klK^) = (46) 

which describes coupled fast, slow and thermal waves. When only longitudinal propagation 
is allowed, the above dispersion relation becomes, 

(cj2 - fc^A^) [^{iriuj + lii) - = (47) 

Then, like in 13.11 we obtain two decoupled Alfven waves, damped by resistivity, whose 
dispersion relation is given by, 

fc' = (48) 
and another dispersion relation, 

= ^ (49) 

which is a fourth degree polynomial in k describing coupled propagating thermal and 
slow waves, only damped by thermal effects. The dispersion relation (j46p has been solved 
numerically and Figures S] and [S] show the behaviour of the damping length, wavelength and 
ratio of the damping length to the wavelength versus period for the fast, slow and thermal 
waves. Starting with Figure|3]we have only considered an interval of periods between 10^^ 
and lO*" s, since for periods smaller than 10^^ s, much shorter than those of interest, the 
curves become very entangled. In the case of the fast wave, when an FIIP is considered the 
spatial damping is governed by radiative losses and thermal conduction (Carbonell et al. 
2006) in the interval of periods from 10~^ to 10*" s. However, in a FIRP we observe a slight 
change in the damping length of fast waves around a period of 1 s. This change tells us 
that the dominance of thermal conduction appears slightly before than in the ideal case. 
Also, looking to the ratio between the damping length and the wavelength we observe that 
in this case the efficiency of the damping for fast waves is improved, with respect to a FIIP, 
for periods below 1 s. For slow waves, no differences appear in the behaviour of FIIP and 
FIRP. The behaviour of thermal waves (Figure [5]) is exactly the same in FIIP and FIRP. 



4.3.2. Partially ionised plasma 

Now, our dispersion relation is given by Eq. (j45l) and when only longitudinal propagation 
is allowed the expression is similar to that of 14.3.11 although the numerical value of the 
Cowling's resistivity is different and the results for fast waves would be different. When 
oblique propagation is considered, the dispersion relation (j45p has been solved numerically 
and a strong distortion of the damping length and wavelength curves corresponding to fast 
waves (Figure 21 left panels) appears. Basically, the changes affect the radiative plateau, 
between periods 10"^ and 10"^ s, and partial ionisation decreases the damping length of 
fast waves in this region. For slow waves (Figure 01 right panels), a similar behaviour is 
found in the same regions although the distortion is not so important since a very short 
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radiative plateau, betvifeen 10^ and 10'^ s remains, together with a region, betvifeen 1 and 
10^ s, where thermal conduction is dominant. As compared to a FIRP, the ratio L^/X 
for fast waves decreases substantially for periods below 10^ s, although for the periods of 
interest in prominences this ratio remains very large. For slow waves, partial ionisation 
produces that the ratio L^/X reaches a maximum of efficiency, ^ 1, for periods similar to 
those involved in prominence oscillations, and this maximum is displaced towards longer 
periods when ionisation is decreased. The changes in the wavelengths of slow and fast waves 
are similar to those shown in the adiabatic case (Sect. HT^ . In the case of thermal waves 
(Figure [5]), partial ionisation increases both the damping length and wavelength of these 
waves, although the behaviour of the ratio id/A is similar to previous cases. Since thermal 
wave is always strongly damped, which makes its detection very difficult, in the following 
we will avoid further comments on it. 

In order to understand the effects of radiation and thermal conduction by neutrals 
and electrons, in Figure IHl we represent, for fast and slow waves, the same quantities but 
with optically thin radiation and heating removed i.e. only thermal conduction is at work, 
and we can observe that the most efficient damping for both waves occurs for partially 
ionised plasmas, which suggests that the inclusion of the isotropic thermal conduction due 
to neutrals plays a very important role for all the periods considered. However, we must 
take into account that when the ionisation fraction decreases, radiation also decreases and 
that, because of neutrals, thermal conduction is favoured, what makes difficult to establish 
meaningful comparisons between plasmas with different degrees of ionisation. 

Next, in Figure [7] we plot the behaviour of fast and slow waves, when thermal conduc- 
tion has been removed i.e. only radiation and heating remain. When thermal conduction 
by electrons and neutrals is removed, the behaviour of fast waves in a FIIP and in FIRP is 
similar and, within the interval of periods of interest, the behaviour of the different quan- 
tities corresponding to the different types of plasma is the same. However, when periods 
below 100 s are considered, FIIP and FIRP are strongly affected by the lack of thermal 
conduction. For slow waves, a similar behaviour appears and the most important conclusion 
is that when we compare Figures [6] and [71 and within the interval of periods of interest, 
the presence of a very efficient damping is due to optically thin radiation. 

On the other hand, when non-adiabatic magnetoacoustic waves are considered, the 
importance of radiation and thermal conduction can be also quantified in terms of two 
dimensionless parameters (De Moortel and Hood 2004), namely, the thermal ratio, which 
modified for the case of partial ionisation becomes. 



which is 1/7 times the ratio of the sound travel time (ts — l/cg) to the thermal conduction 



d = 



(7 - l)(^c|| + Kn)ToPo _ 1 Ts 



(50) 



I'^pIts 7 Tcond 



timescale (rcond 



I'^Po/ [(7 — l)('*ojl + '*n)Tb])i and the radiation ratio 



r = 



(7 - l)rsC^plx*Tl 



(51) 
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which is the ratio of the sound travel time to the radiation timescale (tj- — 
IPo/ [(7 ~ l)^jPoX*2^o'] )■ From the equilibrium parameters, we can compute the value of I 
at which the condition c? = r is satisfied, 

'-^^ 

where Ke\\ and k„ are the thermal conduction coefficients corresponding to anisotropic 
electronic and isotropic, by neutrals, thermal conduction, respectively. Then, when the 
spatial length of the perturbation, the wavelength, is of the order of I or smaller, thermal 
conduction becomes dominant. For the slow wave. Figure |3] (Right panels) shows that the 
transition from a regime dominated by radiation to another dominated by conduction can 
be clearly seen in the plot of the damping length versus period. For a FIIP, this transition 
occurs at a period of 1 s while for a PIP, with fi — 0.8, occurs at a period between 10 and 
100 s, and with fi = 0.95, occurs at a period close to 100 s. Using the values assumed for 
prominence parameters, for FIIP we obtain I 4700 m; for a PIP (// = 0.8), I ~ 40000 m; 
and for a PIP (/i = 0.95), I ~ 100000 m. Then, setting these periods in the plot of the 
wavelength versus period in Figure S] (Right panels), we can check that the numerical value 
of the wavelength is almost coincident with the above analytical determinations for I. The 
reason for this increase in the numerical value of I is that for a PIP thermal conduction is 
enhanced due to neutrals contribution, increasing when the ionisation fraction decreases, 
while the denominator of I decreases when the ionisation fraction decreases. Summarizing, 
when ionisation decreases the period at which the dominant damping mechanism changes 
from radiation to thermal conduction increases and, consequently, the wavelength also 
increases. 



4.4. Non-adiabatic Magnetoacoustic waves in a prominence plasma witti background flow 
4.4.1. Fully ionised ideal plasma 

Setting in Eq. conditions corresponding to FIIP, we obtain, 

{Q^ - lt^A'^)i-Q^) + It'^vKn^ - /fc^A^) = (53) 

which is a sixth degree polynomial in the wavenumber It. When only longitudinal propa- 
gation is allowed, we obtain two undamped Alfven waves given by, 

(54) 

whose solutions for the wavenumbers are 

fc=^^ (55) 

as in Carbonell et al. (2009), while the dispersion relation, 

= k^A^ (56) 

describes coupled slow and thermal waves modified by the flow and damped by thermal 
effects. In the case of oblique propagation, dispersion relation (|53p has been solved nu- 
merically and the behaviour of fast and slow waves is shown in Figure [51 When a flow is 
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present the unfolding of wavelengths and damping lengths appears. Since the considered 
flow speed is much smaller than Alfven speed, the separation of the curves corresponding 
to the fast wave is very small while, since flow speed and sound speed are comparable, the 
curves corresponding to slow waves substantially separate. This effect strongly affects the 
behaviour of the damping length versus wavelength for slow waves, since one of them has a 
very efficient spatial damping for periods observed in prominence oscillations. In the case 
considered in this section, the damping of fast and slow waves is strictly due to thermal 
effects. 

4.4.2. Fully ionised resistive plasma 

In this case, considering FIRP conchtions, the dispersion relation becomes 

{n^ ~ k^K^)(ik^'qVt - r?^) + k^vl(n^ - klK^) = o (57) 

which is a seventh degree polynomial in the wavenumber k. Considering only longitudinal 
propagation, we find, again, coupled slow and thermal waves modified by the flow and 
damped by thermal effects, and three Alfven waves given by the dispersion relation (IMl) . 
with 9 — and its solutions. Therefore, when we solve dispersion relation (j57p . we expect 
three fast, two slow and two thermal propagating waves. In Figure [S] fast and slow waves 
have been plotted and compared with the previous case. Such as it is shown, for fast and 
slow waves no important differences in the behaviour with respect to the ideal case are 
seen. Thermal waves, as well as the third fast wave, are strongly damped after a very short 
distance. 

4.4.3. Partially ionised plasma 

Now, the dispersion relation is given by Eq. (fTO|) and Figure [3] displays the behaviour of 
the damping length, wavelength and ratio of damping length versus wavelength for fast 
and slow waves. The most interesting results are those related with the ratio L^/X. For fast 
waves, this ratio decreases with the period becoming small for periods below 10^^ s while 
for one of the slow waves, the ratio becomes very small for periods typically observed in 
prominence oscillations. When ionisation is decreased, slight changes of the above described 
behaviour occur, the most important being the displacement towards longer periods of the 
peak of most efficient damping corresponding to slow waves. Such as it has been pointed 
out before, the presence of a third fast wave is, again, due to the joint action of flow and 
resistivity. In absence of flow or resistivity, the dispersion relation would become a sixth 
order polynomial in the wavenumber and this wave would be absent. 

5. Conclusions 

Quiescent solar prominences and filaments are partially ionised plasmas and some of their 
typical features are the presence of material flows and oscillations. While the time damp- 
ing of these oscillations has been thoroughly studied, their spatial damping demands an 
in-depth study. Interpreting the observed oscillations in terms of MHD waves, we have anal- 
ysed the spatial damping of Alfven and non-adiabatic magnetoacoustic waves in a flowing 
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partially ionised prominence plasma. Several different cases, with dispersion relations of 
increasing complexity, have been considered, and the conclusions derived from our study 
are summarised in the following. 

As it is well known, Alfven waves are difficult to damp since non-adiabatic effects do not 
affect them, however, when Alfven waves in a partially ionised plasma are considered, they 
can be spatially damped and analytical expressions describing their spatial damping can be 
obtained. When the ionisation decreases, the damping length of these waves also decreases 
and the efficiency of their spatial damping in the range of periods of interest is improved, 
although the most efficient damping is attained for periods below 1 s. A new feature is that 
when a flow is present a new third Alfven wave, strongly attenuated, appears. The presence 
of this wave depends on the joint action of flow and resistivities, since in absence of flow, 
or for a FIIP, the dispersion relation becomes quadratic giving place to the two well-known 
Alfven waves. Furthermore, this third wave could only be detected by an observer not 
moving with the flow. 

When adiabatic magnetoacoustic waves are considered and the effect of partial ionisa- 
tion is taken into account, some new features appear. When a FIRP is considered and only 
longitudinal propagation is allowed, slow waves are decoupled from fast waves, propagating 
undamped while fast waves propagate with a modified Alfven speed and are damped by 
resistivity. When a PIP plasma is studied and only longitudinal propagation is allowed, the 
same happens but now the numerical value of Cowling's resistivity is greater than before 
enhancing the damping. Then, the behaviour of slow waves is only influenced by partial 
ionisation when oblique propagation is allowed. Furthermore, when a PIP is considered the 
behaviour of fast waves is very similar to that of Alfven waves and the damping becomes 
very efficient for periods below 1 s, while for slow waves the peak denoting the most efficient 
damping moves towards higher periods when the plasma ionisation decreases. When in the 
adiabatic case a flow is considered, the main difference is the unfolding of the damping 
length, wavelength and damping length to wavelength curves, and the apparition of a third 
fast wave, strongly damped, due to the join presence of flow and resistivities. 

Considering now non-adiabatic magnetoacoustic waves, when partial ionisation is 
present the behaviour of fast, slow and thermal waves is strongly modified. Comparing 
with non-adiabatic fast waves in a FIIP, which arc damped by electronic thermal conduc- 
tion and radiation, the damping length of a fast wave in a PIP is strongly diminished by 
neutrals thermal conduction for periods between 0.01 and 100 s, and, at the same time, 
the radiative plateau present in FIIP and FIRP disappears. The behaviour of slow waves 
is not so strongly modified as for fast waves, although thermal conduction by neutrals 
also diminishes the damping length for periods below 10 s, and a short radiative plateau 
still remains for periods between 10 and 1000 s. Finally, thermal waves are only slightly 
modified although the effect of partial ionisation is to increase the damping length of these 
waves, just the opposite to what happens with the other waves. Next, when a background 
flow is included, a new third fast wave appears which, again, is due to the joint action of 
flow and resistivities. As we already know, wavelengths and damping lengths arc modified 
by the flow, and since for slow waves sound speed and observed flow speeds are comparable 
this means that the change in wavelength and damping length are important leading to an 
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improvement in the efficiency of the damping. Also, the maximum of efficiency is displaced 
towards long periods when the ionisation decreases, and for ionisation fractions from 0.8 
to 0.95 it is clearly located within the range of periods typically observed in prominence 

oscillations with a value of id/A smaller than 1. This means that for a typical period of 
10^ s, the damping length would be between 10^ and 10"^ km, the wavelength around 10 * 
km and, as a consequence, in a distance smaller than a wavelength the slow wave would be 
strongly attenuated. On the other hand, during our calculations we have found that the 
different heating mechanisms usually considered (Carbonell et al. 2004) do not affect the 
results. 

In conclusion, the joint effect of non-adiabaticity, flows and partial ionisation allows 
slow waves to damp in an efficient way within the interval of periods typically observed 
in prominences. Thermal waves are attenuated very efficiently within the interval of inter- 
est but, probably, their observational detection is very difficult, and fast waves are very 
unefficiently atenuated within the considered interval. Furthermore, since Alfven and fast 
waves display a similar behaviour, the observational discrimination between them should 
be based on the detection of their associated perturbations. It is also important to remark 
that the new fast and Alfven waves are only detectable in a reference system external to 
the flow, however, its short damping length should make its detection very difficult. Such 
as we have seen, even in the most simple case of an imbounded medium threaded by a 
uniform magnetic field, the inclusion of non-adiabatic effects, partial ionisation and flows 
complicates the study of the spatial damping of prominence oscillations because of the 
apparition of new waves and the difficulty to discriminate between the diflFerent effects. 
This also points out that from observations, and due to the entanglement of the different 
effects, it is extremely difficult to properly interpret the observed oscillations in terms of 
MHD waves. 

The present study represents a first step in the investigation of the behaviour of the 
spatial damping of MHD waves in partially ionised prominence plasmas. Forthcoming stud- 
ies must focus in investigating the behaviour of this type of damping in fine structures, 
threads, which we think are the basic constituents of quiescent prominences. 
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